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Abstract 

NGC 4636, an X-ray bright elliptical galaxy, was observed for 70 ks with Suzaku. 
The low background and good energy resolution of the XIS enable us to estimate the 
foreground Galactic emission accurately and hence measure, for the first time, the 
O, Mg, Si and Fe abundances out to a radius of ~28 arcmin (~ 140 kpc). These 
metal abundances are as high as >1 solar within the central 4' and decrease by ~50% 
towards the outer regions. Further, the O to Fe abundance ratio is about 0.60-1.0 
solar in all regions analyzed, indicating that the products of both SNe II and SNe 
la have mixed and diffused to the outer regions of the galaxy. The O and Fe metal 
mass-to-light-ratios (MLR) of NGC 4636 are 2-3 times larger than those of NGC 1399 
implying that metal distributions in NGC 4636 are less extended than those in NGC 
1399, possibly due to environmental factors, such as frequency of galaxy interaction. 
We also found that the MLRs of NGC 4636 at 0.1 rigo are ~5 times smaller than those 
of clusters of galaxies, possibly consistent with the correlation between temperature 
and MLR of other spherically symmetric groups of galaxies. We also confirmed a 
resonant scattering signature in the Fe xw line in the central region, as previously 
reported using the XMM- Newton RGS. 

Key words: galaxies: elliptical and lenticular, cD — X-rays: galaxies — X-rays: 
ISM 
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1. Introduction 



The universe is filled with matter in various physical states. One of these is the hot 
interstellar medium (ISM) in elliptical galaxies. The ISM is gravitationally bound by the mass 
of stars, dark matter and itself, and its temperature is about 0.5-1 keV. Emission from the ISM 
is observed in the X-ray band and its luminosity reaches about 10 39 -10 42 erg s -1 . By accurate 
analysis of X-ray spectra, we can estimate the temperature and metallicity distributions of the 
ISM and reveal the metal enrichment history of the universe. 

In the past, X-ray imaging detectors, such as Chandra ACIS or XMM-Newton EPIC, 
have observed the ISM with high angular resolution and high signal-to- noise ratio (S/N) or wide 
field of view, and revealed temperature structures, surface brightness and mass distributions 
(e.g., Buote 2002; Buote et al. 2003; Fukazawa et al. 2006; Humphrey et al. 2006). In addition, 
observations by XMM-Newton have suggested similar metal abundance patterns to those in the 
centers of clusters or groups of galaxies (e.g., Matsushita et al. 2003; Matsushita et al. 2007b). 
The excellent energy resolution of XMM-Newton RGS has been used to precisely measure metal 
abundances (e.g., Xu et al. 2002; Tamura et al. 2003) and show that Si and Fe abundances 
in the central region are similar, while the O abundance is in general about half that of Fe. 
To understand metal enrichment processes, measurements of O, Ne, or Mg abundances in the 
outer regions are needed. To unambiguously constrain the O abundance in elliptical galaxies it 
is also necessary to have an accurate estimate of the foreground Galactic thermal emission. The 
best detector to satisfy these conditions is the Suzaku XIS which has a better energy resolution 
and a lower background than any previous X-ray CCD camera. The launch of the Suzaku XIS 
(Koyama et al. 2007) with high energy resolution in the lower energy X-ray band allows us 
to obtain O and Mg abundance distributions through accurate measurements of their emission 
lines. Past studies (e.g., NGC1399 ; Matsushita et al. 2007a, NGC 5044 ; Komiyama et al. 
2008) have shown that the O to Fe abundance ratio, which can be used to estimate the number 
ratio of la to II supernovae, is about 0.6-0.8 solar in all regions. However, observations in the 
outer regions, where metals ejected from galaxies have remained, are still poor. So far, Suzaku 
measurements of abundance distribution have been reported for three elliptical galaxies. NGC 
5044 and NGC 507 exhibit a steeper radial O abundance profile than NGC 1399 while there 
is a large scatter among their O mass-to-light ratios. (Matsushita et al. 2007a; Komiyama et 
al. 2009; Sato et al. 2009a). Studies of these properties in additional galaxies will thus aid 
understanding of metal enrichment processes. 

In this paper, we report the results of spectral analysis of NGC 4636 observed by the XIS 
detector onboard Suzaku (Mitsuda et al. 2007). NGC 4636 (z = 0.00313; ~ 17 Mpc, Smith et al. 
2000) is an X-ray bright elliptical galaxy belonging to the Virgo cluster. Metal-enriched hot gas 
is extended over the entire galaxy group (Matsushita et al. 1998). This galaxy is located at the 
south end of the cluster, where the X-ray emission from the Virgo center is negligible, so we can 
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study the metal enrichment history of a giant elliptical galaxy (Matsushita et al. 1997). So far, 
observations by ASCA (Matsushita et al. 1997), Chandra (Ohto et al. 2003; Humphrey et al. 
2006), and XMM- Newton (e.g. Xu et al. 2002) have been used to measure metal abundances. 
The temperature of hot gas in this galaxy is 0.6-0.7 keV, which is lower than the 1.0-1.5 
keV of NGC 1399, NGC 507, and NGC 5044. Therefore, we can also discuss the temperature 
dependence of metal properties. In addition, NGC 4636 is quite isolated in comparison with 
the above 3 galaxies, and thus the comparison will aid investigations of environmental effects. 
We describe the observation and data reduction in section 2. Section 3 on the analysis and 
results describes the methods we used to estimate background and fit spectra as well as the 
resulting temperature and metallicity distributions. Additionally, we discuss the confirmation 
of resonance scattering in the Fe xvn line in the spectrum of the central region. In section 4, 
we discuss the O, Ne, Mg, Si abundance ratios to Fe, their distributions, the metal transport 
process, and resonance scattering of the Fe xvn line. A summary is given in section 5. 

Throughout this paper, we use a Hubble constant value of H = 70 km s _1 Mpc _1 or /t 10 o 
= 0.7 (1 arcmin ~ 5 kpc), and a virial radius of ri 8 o = 1.95/i^jo \fk(T) /10keV Mpc (Markevitch 
et al. 1998) where kT is the gas temperature. The Galactic hydrogen column density is A% 
= 1.8xl0 20 cm" 2 in the direction of NGC 4636 (Stark et al. 1992). We use the Anders and 
Grevesse (1989) definition of Solar abundance ratios for easier comparison with other Suzaku 
papers which use this definition. Other recent papers often use instead Grevesse and Sauval 
(1998). The former has 26% and 48% higher abundances for O and Fe, respectively. All errors 
are given as a 90% confidence range. 

2. Observation and Data Reduction 

A pointing observation of NGC 4636 by Suzaku was carried out on December 6-7, 2005 
for about 70 ks. In order to measure metal abundances in the outer region of NGC 4636, a 
10' north offset observation was carried out on December 7-9, 2007. Another north-west offset 
observation at 3.2° away from the NGC 4636 center was also carried out on June 17-18, 2007 to 
estimate the foreground Galactic emission - this region contains no other bright objects. All the 
XIS data were taken with both 3x3 and 5x5 pixel mode. The observation log is summarized 
in Table 1. 

We analyzed the public rev 1.2 data and screened them with cut-off rigidity (COR) > 8 
GV, target elevation angle above the earth rim (ELV) > 5°, and elevation angle above the day 
earth rim (DYE_ELV) > 20°. For the 2005 observation, we extracted spectra of six annular 
regions centered on the X-ray peak (Figure 1(a)). For the north offset observation, we analyzed 
the spectrum of a 8'-radius aperture centered on a position at 20' away from the galaxy center 
(Figure 1(b)). In extracting the above spectra, we excluded the region which is contaminated by 
the Fe calibration source. We used xisrmfgen version 2006-11-26 and xissimarfgen version 2006- 
10-26 to make response matrix files (RMF) and ancillary response files (ARF), respectively. The 
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Table 1. Observation log 



rarget 


iNlcrO 4DOD 


lArrr 1 ac\ic\ \i o iy r i i 

40o0-1n Uivl rl 




z 


0.00313 






position (Ra, Dec) in J2000 


190.7250, 2.7520 


190.6828, 3.0503 


192.5138, 5.4608 


exposure time (ks) 


70.3 


51.1 


34.9 


date (y-m-d h:m) 


05-12-06 06:38 


07-12-07 05:01 


07-06-17 20:11 



^GC 4636-NORTH : 10' offset observation to the north of NGC 4636 

2 NGC 4636-GALACTIC : 3.2° offset observation for estimating of the foreground Galactic emission 



source mode of xissimarfgen was set to SKYFITS, which is appropriate for analyzing extended 
objects, and the XIS image file was input. Data from the hard X-ray detector, HXD, are not 
treated here because no signal was detected from NGC 4636; an upper limit of the hard X-ray 
emission is 9 x 10~ 12 erg cm" 2 s" 1 (10-50 keV) for the 3% background uncertainty (Fukazawa 
et al. 2009). 



(a) (b) 




Ra(J2000) Ra(J2000) 



Fig. 1. X-ray image (0.3-5keV) of NGC 4636, convolved with a gaussian of a = 3 pixel. The NXB 
(non X-ray background) is subtracted, and the vignetting effects and the difference in exposure times arc 
corrected, (a) Observation of the center region, (b) Mosaic image of three observations. Boundaries of 
analysis regions are shown by a circle. 



3. Spectral Analysis and Results 

3.1. Estimation of Background 

In order to study the emission from NGC 4636 correctly, we have to estimate the back- 
ground accurately. This consists of the instrumental non X-ray background (NXB), the cosmic 
X-ray background (CXB) and the foreground Galactic X-ray emission (GXE). 

First, we estimated the NXB spectra in each analysis region using the tool "xisnte- 
bgdgen", which adds COR-sorted (< 4GV, 4-13 GV in 1 GV step, and > 13 GV) night-Earth 



4 



Table 2. Spectral parameters of the CXB and the GXE 



Cosmic X-ray Background (CXB)* 


POWER-LAW x HIGHECUT photon-index 


normalization^ 


e-folding energy (keV) 


1.29 


8.21xlCT 4 


40.0 


Galactic X-ray Emission (GXE)t 


APEC + APEC Milky Way Hallo (MWH) 


Local Hot Bubble (LHB) 




kT (keV) normalization^ 


kT (keV) normalization^ 


reduced-^ 2 (d.o.f.) 


0.326±0.013 2.31±0.37 


0.094±0.004 2.02±0.84 


1.48 (206) 



* Parameters of the CXB are fixed to the values obtained with HEAO-1 (Boldt ct al. 1987). 
f Parameters of the GXE are from the best fit result for spectra of the 3.2° offset region. 



Metal abundance of the APEC model is fixed to 1 solar. 
\ normalization = photons cm -2 s _1 keV -1 (7r20 2 arcmin)" 2 QlkcV. This normalization is estimated 

by assuming in the uniform-sky ARF calculation (20' radius). 
§ normalization is for the APEC model and equivalent to 1{)~ 17 / {An(l+z) 2 D 2 A } j n c iiH^cm" 5 

(Da '■ angular size distance to the source). 

spectra weighted by the exposure time for each COR in the NGC 4636 observation. Second, 
using the sky-averaged CXB parameters from HEAO-1 (Boldt et al. 1987) summarized in table 
2, we estimated the CXB contribution by using an uniform-sky ARF, which is generated with 
"xissimarfgen" and assumes that the emission is uniform within 20'. Finally, we estimated the 
GXE, using the 3.2° offset observation by extracting spectra within 8' (~ 40 kpc) of the center 
of the field of view. We used a uniform-sky ARF as for the CXB estimation. The GXE consists 
of a Local Hot Bubble (LHB) component and a Milky Way Halo (MWH) component. These 
have temperatures of approximately ~0.08 keV and ~0.2-0.3 keV, respectively (e.g., Lumb et 
al. 2002). We fitted the GXE spectrum with a two-temperature APEC model, where temper- 
atures are set to be free and metal abundances are fixed to one solar. At these temperatures, 
the continuum emission is vanishingly small compared to the line emission. Therefore, even 
if this assumption for metal abundance is not completely correct, the estimation of the GXE 
will not change significantly. The results of the fit are shown in table 2, along with the fixed 
parameter values used for the CXB model. Figure 2 shows the best fit model and spectra. The 
temperatures obtained are consistent within errors with the results from the ROSAT-PSPC all- 
sky survey (Kuntz & Snowden 2000). Using the best-fit parameters, we modeled the Galactic 
emission in each spectral region of the NGC 4636 observations by assuming the same surface 
brightness and spectral parameters as those of the 3.2° offset region. 

3.2. Spectral Fitting 

We fitted the spectra of each region of NGC 4636, based on the above estimation of 
the NXB, the CXB, and the GXE. The CXB and the GXE are included as models in the 
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Fig. 2. Best-fit result of spectral analysis of the Galactic emission for the 3.2° offset observation. Black 
and red crosses are observed spectra of FI and BI, respectively. Black and red thick lines are the best-fit 
models of FI and BI, respectively. As for individual model components, solid lines are for FI and dot lines 
are for BI. Orange is the CXB model, purple and sky-blue are the Galactic emission model for MWH and 
LHB, respectively. Residuals are shown in the bottom panels. 

spectral fit with their spectral parameters fixed, and the NXB was subtracted. We applied a 
two-component model for the emission from NGC 4636 : a single-temperature vAPEC model 
for the ISM, and a BREMSS model with a temperature fixed to 7 keV (Matsushita et al. 1994) 
for the hard component due to low- mass X-ray binaries (LMXB). For the all regions, the fitting 
is performed in the energy range of 0.3-5 keV because the spectrum below 0.3 keV is difficult to 
analyze due to large systematic errors, and the NXB dominates above 5 keV. The energy range 
around the Si K-edge (1.82-1.841 keV) was ignored due to response problems. We divided the 
metals into several elemental groups (He = C = N, O, Ne, Mg = Al, Si, S = Ar = Ca, Fe 
= Ni), and fixed the abundance of the He=C=N group to solar. Abundances of other groups 
were treated as free parameters. 

Table 3 shows the best-fit results for each region. The fit is acceptable for all the regions. 
The O abundance is measured up to 10-12' and Mg, Si and Fe to 12-28' for the first time. 
However, the Ne abundance is not reliable because Ne line measurements are influenced by 
nearby Fe-L lines, which cannot be resolved with the resolution of the XIS. It can be seen 
from Figure 3 that the Galactic emission is well reproduced and significantly contributes to the 
spectra in the outer region, and thus it must be accurately estimated to derive the O abundance 
of NGC 4636. As a result, even for O, the abundances in the outer regions are determined up to 
10' with ~30% accuracy. Even if we allow the temperature of the BREMSS (7keV) component 
to vary, the metal abundances only change within errors so we believe that they are reliable. 
The luminosity of the hard X-ray component of the BREMSS model is 1.5 xlO 40 erg s _1 over 
2-10 keV within 6' of NGC 4636. This is consistent with the ASCA results (Matsushita et 
al. 1994), indicating that the hard X-ray component is due to LMXB. We also confirmed that 
the abundance results are almost the same when we fit the hard component with a power-law 
of index 1.5, which was used for the LMXB estimation in Loewenstein et al. (2001). In the 
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Fig. 3. Spectral fitting of each annular region of NGC 4636. Black and red crosses are observed spectra 
of FI and BI, respectively. Black and red thick lines are the best-fit models to the FI and BI, respectively. 
For individual model components, solid lines are for FI and dot lines arc for BI. Green and blue are the 
thermal component and LMXB component from NGC 4636, respectively. Orange is the CXB model, 
purple and sky-blue are the Galactic emission models for the MWH and LHB, respectively. Residuals are 
shown in the bottom panels. Shown are spectra of 0-2', 2-4', 4-6' regions in the upper panels, 6-8', 8-10', 
10-12' in the middle panels, and 12-28' in the bottom. 

central region, the reduced-;^ 2 is a little larger than that in the outer region. This might be 
due to temperature structure in the center region (Ohto et al. 2003) so we tried fitting with 
a two-temperature model only in the 0-2' region. We assumed that the metal abundances 
are the same for both temperature components. Comparing the result with that of the single- 
temperature model, the metal abundances do not greatly change and the temperature of the hot 
component in the two-temperature model is almost the same as that in the single-temperature. 
The F-test probability is 0.097% and thus the improvement is marginal. Therefore, hereafter 
we focus on the results with the single-temperature model. A fit residual around 0.7 keV in 
the 0-2' region is discussed and analyzed in §3.4. 

3.3. Temperature and Metal Abundances 

Distributions of temperature and metal abundances are plotted against radius (kpc) in 
Figures 4 and 5, respectively. Systematic errors due to the GXE, CXB, and NXB are also 
shown. The systematic errors due to the GXE are estimated by varying the normalizations of 
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Table 3. Results of spectral fitting of each annulus region of NGC 4636 



region 


model kT (keV) 


O (solar) 


Ne (solar) 


Mg (solar) 


Si (solar) 


Fe (solar) 


0-2' 




0.642±0.002 


0.75±0.06 


1.54±0.12 


1.50±0.11 


1.51±0.12 


1.14±0.07 


0-2' 


2T b 


1.56±0.33, 0.641±0.002 0.79±0.07 


1.47±0.15 


1.52±0.11 


1.54±0.12 


1.20±0.07 


2-4' 


IT 


0.733±0.004 


1.01±0.13 


2.14±0.28 


2.33±0.24 


1.81±0.19 


1.41±0.12 


4-6' 


IT 


0.792±0.005 


0.41±0.06 


0.62±0.11 


0.86±0.08 


0.57±0.06 


0.49±0.03 


6-8' 


IT 


0.811±0.006 


0.41±0.08 


0.45±0.13 


0.74±0.10 


0.50±0.07 


0.43±0.03 


8-10' 


IT 


0.796±0.008 


0.34±0.10 


0.60±0.17 


0.86±0.14 


0.53±0.11 


0.43±0.05 


10-12' 


IT 


0.774±0.013 


0.58±0.28 


0.67± 0.42 


1.20±0.40 


0.84±0.39 


0.68±0.17 


12-28' 


IT 


0.746±0.016 


<1.20 


<1.60 


1.13±0.41 


0.87±0.31 


0.63±0.17 


region 


c vAPEC normalization 


d BREMSS normalization 


reduced- 


X 2 (d.o.f.) 


0-2' 




2.12±0.21 




9.58±1.01 




1.40 (829) 


0-2' 




0.10±0.06, 2.01±0.20 




7.73±2.13 




1.38 (827) 


2-4' 




1.08±0.17 




9.34±1.30 




1.29 (844) 


4-6' 




1.51±0.14 




3.98±1.17 




1.11 (518) 


6-8' 




0.97±0.13 




4.04±1.12 




1.08 (446) 


8-10' 




0.72±0.13 




5.34±1.19 




1.01 (371) 


10-12' 




0.27±0.13 




3.18±1.21 




0.96 (210) 


12-28' 




0.65±0.33 




3.60±2.36 




1.01 (586) 



a single-temperature vAPEC model 

b two-temperature vAPEC model 

c [10- 17 /{47r(l+z) 2 Z)i}]/n c nHrfycm- 5 

(rie : electron density, nu : hydrogen density, Da '■ angular size distance to the source) 
n e nidV cm (m : ion density, D : distance to the source) 



the APEC model for the MWH and LHB within their error ranges, which are comparable to 
the X-ray count fluctuation (15%) of the ROSAT All-sky X-ray background survey (band 4) 1 
for the NGC 4636 region. The systematic errors due to the CXB and NXB are at most 5% 
(Tawara et al. 2008), respectively, and they are dominant in the hard X-ray band. Therefore, 
the systematic errors on the temperature and metal abundances are mostly caused by the GXE 
uncertainty. A 15% uncertainty in the estimate of the XIS optical blocking filter thickness 
introduces systematic errors comparable to those from the GXE. The temperature in the cen- 
tral region is 0.64 keV, rising to 0.7-0.8 keV. This profile is consistent with that from ASCA 
(Matsushita et al. 1997), XMM-Newton (Xu et al. 2002) and Chandra (Ohto et al. 2003) and 
the temperature distribution is barely affected by any systematic errors in the GXE. For the 
inner 50 kpc most metal abundances are tightly restricted thanks to the high energy resolution 
of the Suzaku XIS while in the outer regions the abundances are less well constrained due to 
the GXE systematic errors. However, the Ne abundance is not reliable because Ne lines are 

1 SkyView: http://skyview.gsfc.nasa.gov/ 
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confused with Fe-L lines. Furthermore, only upper limits could be obtained on the O and Ne 
abundances in the 60-140 kpc region, since their uncertainties are very large due to contam- 
ination by the Galactic emission. It can been seen that all the metal abundances are 0.7-2.5 
solar in the central region, and decrease down to 0.3-1.0 solar towards the outer region. In 
the central region, the O and Fe abundances are about 1.2 times larger than that obtained 
with the RGS (Xu et al. 2002). Figure 6 shows confidence contours of O and Mg abundances 
against the Fe abundance. The O and Mg to Fe abundance ratio are particularly important to 
constrain the contribution of different types of supernovae. Although the individual errors for 
each metal abundance are large, the abundance ratios are constrained well, especially for the 
Mg to Fe abundance ratio in the outermost region. 
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Fig. 4. Temperature distribution of NGC 4636 (red) plotted against the radius in units of kpc. Dotted 
lines indicate the systematic error range due to the GXE uncertainty. For comparison, Suzaku results of 
NGC 1399 (blue) (Matsushita et al. 2007a) are also shown. 



3.4. Resonance Scattering in the Central Region 

In Figure 7, we show the spectra in the 0-2'and 2-4' regions. In the spectrum of the 0-2' 
region, the fit residual has a negative line feature around 0.85 keV. This feature disappears in 
the 2-4' region and other outer regions (Figure 3). This is possibly due to resonance scattering 
of the Fe xw line, which was reported with the RGS (Xu et al. 2002; Werner et al. 2009). When 
we introduce a negative gaussian line in the spectral model, the fit significantly improves with a 
F-test probability of < 0.01%. In that case, the Fe abundance increases by 15%. We compared 
the XIS result in the central region with that using the RGS and looked into whether the 
features are consistent. We modeled the XIS spectrum with two BREMSS for the continuum 
thermal emission from the plasma and the LMXBs, and add each emission line with a gaussian. 
The temperatures of these BREMSS models are fixed to 0.642 keV and 7.0 keV for the thermal 
and LMXB emission, respectively. The former is the best-fit vAPEC temperature in table 3. 
The lines considered here are listed in Table 4. Since the XIS cannot resolve individual lines, 
we summed the flux of neighboring lines within 30 eV in the RGS paper (Xu et al. 2002), and 
assumed the strength-weighted line energy. Furthermore, we refer to the ATOMDB line table 
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Fig. 5. Metal abundance distributions of NGC 4636 (red) plotted against the radius in units of kpc. 
Dotted lines indicate the systematic error range due to the GXE uncertainty. For comparison, Suzaku 
results of NGC 1399 (blue) (Matsushita et al. 2007a) are also shown. The O and Ne abundances in 60-140 
kpc represent an upper limit. 

2 for the energy and intensity of the 0.961 and 0.963 keV lines, which were not analyzed in the 
RGS spectrum, but are significant in the XIS spectra. The Ne lines around this energy range 
are not shown because their fluxes are not reliable due to the influence of Fe-L lines. The line 
width of all lines is fixed to keV. 

The line intensity ratio against the summed flux of the strongest line (0.737 and 0.728 
keV ) of Fe xvn is summarized in Table 4 and Figure 8. In Figure 8, the predicted line intensity 
ratios are also shown. These were obtained by simulating spectra with the emission models 
and parameters from Table 3 and fitting them with the same model used for the data. The 
result indicates that the flux intensity of the Fe xvn line (0.825 keV), which was suggested 
to be affected by resonance scattering (Xu et al. 2002), is about 0.75 times lower than that 
predicted by the model. This might indicate that the resonance scattering occurs for the Fe 
xvn ion. The Suzaku intensity ratio is about 1.4 times higher than that from the RGS. This 
might be because the RGS photons are accumulated from a smaller central region than the 
XIS where the resonance scattering is more effective. Doron and Behar (2002) indicated the 
problem of modeling the Fe+16/Fe+17 ion lines, but their recalculation showed that the line 
intensity ratio of the 0.825 keV line to the sum of the 0.737 and 0.728 keV lines changes by 
30-40%, mostly independent of the temperature in the range of 0.6-0.9 keV. Their paper does 
not refer to the APEC model, and thus this factor of change is not necessarily accepted for 
our cases. Our results show that the discrepancy of the line intensity ratio between the data 



2 http : / / cxc . harvard . edu/atomdb/WebGUIDE/ index . html 
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Fe (solar) Fe (solar) 



Fig. 6. Confidence contours of Fe abundance vs O abundance (left), and Fe abundance vs Mg abundance 
(right) for the inner (top) and outer (bottom) region. The square mark denotes the best-fit location, and 
the contours represent 68%, 90%, 99% confidence level, respectively. 

and the model is seen only at the central region. Therefore, we conclude that the resonance 
scattering is robust. 




channel energy (keV) channel energy (keV) 

Fig. 7. Spectrum of 0-2'(left) and 2-4'(right) region. Black and red crosses are observational spectra of 
FI and BI, respectively. Black and red lines are the best-fit models of FI and BI, respectively. Residuals 
are shown in the bottom panels. 
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Table 4. Summary of line intensities in 0.4—1.5 keV for the RGS (Xu et al. 2002) and the XIS. Line intensity ratio is against 
the sum flux of 0.737 keV and 0.728 keV Fe lines. 



Ion 






RGS 






XIS 






E (keV) 


*fiux 


*sum of flux 


intensity ratio 


E (keV) 


*flux 


intensity ratio 


N vi 


0.418-0.428 


0.43 


0.43 


0.05 


0.428 








N vn 


0.497 


0.20 


0.20 


0.03 


0.497 








vn 


0.557 


0.21 












vn 


0.565 


0.09 












vn 


0.570 


0.13 


0.43 


0.05 


0.563 








vm 


0.651 


1.96 


1.96 


0.245 


0.651 


1.03±0.18 


0.21±0.04 


Fe xvn 


0.728 


5.66 












Fe xvn 


0.737 


2.33 


7.99 


1 


0.731 


4.85±0.26 


1.00±0.08 


Fe xvm Fe xdc vm 


0.769 


1.99 


1.99 


0.249 


0.769 


2.94±0.28 


0.61±0.07 


Fe xvn 


0.809 


1.66 












Fe xvn 


0.825 


3.56 


5.22 


0.653 


0.820 


4.46±0.26 


0.92±0.07 


Fe xvm 


0.853 


2.25 












Fe xvm 


0.859 


1.85 












Fe xvm Fe xx 


0.865 


1.45 












Fe xvm 


0.871 


1.90 


7.45 


0.932 


0.861 


4.98±0.23 


1.03±0.07 


Fe xx 


0.961 














Fe xx 


0.963 








0.962 


1.73±0.15 


0.36±0.04 


Mg xi 


1.35 


0.59 


0.59 


0.070 


1.35 


0.61±0.04 


0.13±0.01 



*fiux : 10 4 pliotons cm 2 s 



4. Discussion 

We measured the temperature and metal abundance distributions of NGC 4636 out to 
28' and confirmed the resonance scattering in the central region. In this section we compare 
the Suzaku results with those for other X-ray bright elliptical galaxies (NGC 1399, NGC 507, 
and NGC 5044), and discuss the metal enrichment process in elliptical galaxies. 

4-1. Metallicity Distribution 

In Figures 4 and 5, observational results by the Suzaku XIS for NGC 1399 (Matsushita 
et al. 2007a) are also shown. These are the results with a two-temperature model within 4' 
and one-temperature model beyond 4'. The temperature of NGC 1399 is about twice as high 
as that of NGC 4636, indicating that the gravitational potential of NGC 1399 is deeper. In 
Figure 5, the Ne abundance of NGC 1399 was not measured because its line was polluted by 
Fe-L lines. The radial distributions of the O, Mg, Si and Fe abundances are similar between 
the two galaxies, but metal abundances are about twice as high as those of NGC 1399 in all 
radial regions. In Figure 9, radial profiles of the abundance ratios of O, Ne, Mg, and Si against 
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Fig. 8. Relative line intensity against the sum flux of 0.733 keV and 0.728 keV lines of Fe xvn- Squares 
and circles are the data of the XIS and RGS (Xu et al. 2003). The horizontal axis represents line energy, 
and the vertical axis relative line intensity. Model prediction (filled triangle) is obtained by simulating 
spectra with the emission models and parameters on Table 3 and fitting them with BREMSS and ZGAUSS 
models. 

Fe are shown, with error bars estimated from the allowed slope of the confidence contours in 
Figure 6. For comparison, the abundance ratios for a model of SN II (for all elements) and 
SN la (for Si) are also shown. The abundance ratios for SN II are calculated by integrating 
the metal mass ejected by SNe II (Iwamoto et al. 2006) with a progenitor mass from 1OM 
to 50M Q , weighting by the Salpeter initial mass function, and calculated, for example, as O: 
1.80M© and Fe: O.O9O7M . The W7 model (O: 0.143M o and Fe: O.749M ; Iwamoto et al. 
1999) is used for SN la. The result indicates that all metal abundances can be explained by a 
mixture of SN II and SN la products. This indicates that, in NGC 4636 and NGC 1399, each 
metal diffused to the outside of the galaxy after mixing products from both SN II and SN la. 
For NGC 4636, all the abundance ratios are constant in the radial direction. The abundance 
ratios against Fe in the central region (0-10 kpc) are 0.66±[};Qg solar, 1.35±q;}q solar, 1.32±q;q7 
solar, 1.32d=Q ^Q solar for O, Ne, Mg, Si, respectively. The O/Fe ratio is similar to that of NGC 
1404 (Matsushita et al. 2007a) and NGC 720 (Tawara et al. 2006) measured by the XIS, and 
also with the RGS results for NGC 4636 (Xu et al. 2002) and NGC 5044 (Tamura et al. 2003). 
This indicates that a common metal production process has taken place in elliptical galaxies. 
However, the Mg/Fe ratio of NGC 4636 is higher than that of NGC 1399 in the outer region. 
This difference might be caused by different populations of SN II progenitor stars between the 
two galaxies, but problems with the plasma emission model cannot be ruled out (Matsushita 
et al. 2003). 

It is noticeable that the radial profile of the O abundance of NGC 4636 exhibits a steeper 
gradient than other elliptical galaxies, NGC 1399 , NGC 5044, and NGC 507. This indicates 
that NGC 4636 keeps metals more tightly concentrated. NGC 4636 is rather isolated and has 
never been identified as a galaxy group; the other galaxies are classified as the center galaxy 
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of their galaxy group, which contains many dwarf galaxies (Ferguson and Sandage 1990; Tifft 
et al. 1975). Therefore, the differences in the abundance gradient could be explained by NGC 
4636 experiencing fewer galaxy interactions leading to reduced diffusion of metals to the outer 
region. We consider this further in the next section. 
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Fig. 9. Radial profiles of metals to Fe abundance ratios for NGC 4636 (red) and NGC 1399 (blue; 
Matsushita et al. 2007a). The dashed line represents the ratios for the SN II model (Nomoto et al. 
2006), and the dotted line the SN la model (Iwamoto et al. 1999). The O, Ne and Mg abundance ratios 
for SN la are not shown because these values are quite small (~ 0.01). O and Ne abundances in 60-140 
kpc are upper limits. 



4-2. Metal-mass-to-light-ratio 

The gas mass distributions of NGC 4636 and NGC 1399 are shown in Figure 10, which 
also shows the Chandra results within 0.022 r 180 (— 15.6 kpc) (Fukazawa et al. 2006) and the 
ASCA results beyond this radius (Matsushita et al. 1998). These results show that the gas 
mass of NGC 4636 is about half as much as that of NGC 1399. This is possibly because NGC 
1399 sits at the center of the Fornax cluster, where the potential is deeper than that for NGC 
4636 (Figure 4), so the gas is more tightly bound by NGC 1399. 

The NGC 4636 metal mass distributions and metal-mass-to-light-ratios (MLR) of O and 
Fe are shown in Figures 11(a) and (b), respectively. The metal mass is obtained by multiplying 
the metal abundance value in the outermost region in Figure 5 with the gas mass distribution 
(Figure 10). The stellar light distribution is quoted from Fukazawa et al. (2006). 
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Fig. 10. Gas mass distributions of NGC 4636 (red) and NGC 1399 (blue). For NGC 4636, data within 
and beyond 0.022 ri§o (— 15.6 kpc) refer to the results of Chandra (Fukazawa et al. 2006) and ASCA 
(Matsushita et al. 1998), respectively. For NGC 1399, data within and beyond 0.014 r 180 (~ 14.7 kpc) 
refer to the results of Chandra (Fukazawa et al. 2006) and ASCA (Ikebe 1995) respectively. Double-/? 
model is used as a gas density distribution. 
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Fig. 11. (a) Mass distributions of O and Fc in the hot gas of NGC 4636 and NGC 1399. (b) 
Metal-mass-to-light-ratios of NGC 4636 and NGC 1399. The stellar light distribution is quoted from 
Fukazawa et al. (2006). OMLR and IMLR represents oxygen-mass-to- light ratio and iron-mass-to-ligh 
ratio, respectively. 

Figure 11(a) shows that the O and Fe mass distributions are similar between the two 
galaxies and the metal mass of NGC 4636 is almost the same as that of NGC 1399 at all radii. 
Combining this with the difference of the gas mass, it can be understood that the abundance 
and the gas mass of NGC 4636 are about twice and half those of NGC 1399, respectively. As 
seen in Figure 11(b), MLR of NGC 4636 is also 2~3 times larger than that of NGC 1399 in 
all regions for both metals. This indicates that the metals in the ISM are distributed more 
compactly in NGC 4636; the metals in NGC 4636 do not diffuse so widely. This difference 
might be caused by environmental effects; galaxy interaction is frequent for NGC 1399 and 
thus metal-rich gas could be stirred by galaxy motions. 

Finally, we compared the MLR of NGC 4636 with Suzaku results of other groups and 
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clusters. These objects are galaxy groups centered on NGC 1399, NGC 5044, and NGC 507. 
They are X-ray bright elliptical galaxies like NGC 4636, but the galaxy density is high around 
them and thus they are identified as centers of galaxy groups. We also compare with galaxy 
clusters such as Abell 262, Abell 1060, AWM 7, and Centaurus. Following Sato at al. (2009b), 
we plot the correlation between temperature and MLR at 0.1 rigo in figure 12. This result 
shows that MLR of groups and galaxies with a temperature of ~1 keV, including NGC 4636, is 
~5 times lower than that for clusters. Among galaxies and groups, NGC 4636, NGC 5044 and 
HCG 62 exhibit a spherically symmetric X-ray emission, and NGC 4636 has the lowest MLR 
among them. A possible correlation between temperature and MLR is suggested only for these 
spherically symmetric objects. Such a correlation could be explained by higher temperature 
systems possessing a deeper gravitational potential which prevents dispersal of the metals. 
Although galaxy velocity dispersion is larger for higher temperature clusters, there are no big 
galaxies around cD galaxies. Therefore, we think that the deep potential can confine the metals 
even if a small galaxy passes around the cD galaxy. On the other hand, the X-ray asymmetric 
groups NGC 1399 and NGC 507 show a relatively low MLR. These imply that NGC 1399 
and NGC 507 experience galaxy interactions frequently and metal-rich gas could be stirred by 
galaxy motions. 
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Fig. 12. IMLR and OMLR at 0.1 ri 8 o using B band luminosity. References : NGC 5044 (Komiyama 
ct al. 2008), NGC 1399 (Matsushita ct al. 2007a), NGC 507 (Sato ct al. 2009a), HCG 62 (Tokoi et al. 
2008), Abell 262 (Sato ct al. 2009b), Abell 1060 (Sato ct al. 2007a), AWM 7 (Sato ct al. 2008), Centaurus 
(Matsushita et al. 2007b) 



4-3. Optical depth of Fe xvn the central region 

We calculate the optical depth of resonance scattering for the Fe xvn hne in the central 
region mentioned in §3.4. We assume that the gas density distribution follows the (3 model. 
The scattering cross section a v is represented by the following formula (Shigeyama et al. 1998), 
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,2 



ne mi 

°v = gf\ ttt^ 1 

m e z/ V l-KKl 

rrii is the mass of the ion, gf is oscillator strength, z/ is the frequency of the emission and T 
is the temperature of the system. We set T to 0.642 keV which is the temperature in the 0-2' 
region and the value of gf to 2.5 (Mewe et al. 1985). As a result, a v becomes 1.20xl0~ 15 cm 2 . 
For the density distribution n(r), we use the /3 model, 

n(r) = n [l + (-) 2 ]- |/3 (2) 

77-o is the core density and r c the core radius. For NGC 4636, the core density, the core radius 
and j3 are 0.2 cm -3 , 1.1 kpc and 0.47, respectively (Ohto et al. 2003). 

Then, the optical depth for the emission line from the circular region with projected 
radius R is calculated by the following formula, assuming that the Fe abundance in the central 
region Ap G is 1.14 solar from the value summarized in Table 3, the Fe abundance solar ratio 
Fe/H is 4.68 xlO" 5 (Anders & Grevesse 1989), ionic fraction f\ of Fe xw is 0.23 (SPEX) and 
gas exists up to 300 kpc. 

r(r) = / A Fe (Fe/H)/^ n(0 i dr (3) 

JR J r 2 _ R 2 



As a result of this calculation, the value of r in the central region within 2' becomes 24.1 < r 
< 56.3. Note that those values are not influenced by thermal Doppler broadening and extent 
of the diffuse source. Even if we consider them, r 3> 1 is secure. Our result is somewhat larger 
than that of Werner et al. (2009), and the difference is possibly explained by the assumption of 
the Fe abundance. Therefore, it is optical thick there and resonance scattering will easily occur. 
Turbulence could reduce the resonance scattering (Werner et al. 2009), but the scattering is 
significant for NGC 4636. 

5. Summary 

In this paper, we performed spectral analysis out to 28'(~ 140 kpc) from the center of 
NGC 4636, using data from the XIS detectors onboard Suzaku. We obtained the following 
results. 

• The temperature distribution shows a smooth gradient from the center (0.642 keV) to the 
outer regions. For the metal abundance, we found that the distributions of O, Mg, Si, Fe 
are all higher in the central 4' region and decrease by 50% toward the outer region. It is 
noticeable that O and Mg abundances are determined accurately, even out to 28' from the 
center. 

• The abundance ratio of O/Fe is about 0.6~1.0 solar independent of the radius. This means 
that O and Fe in the gas diffused to the outside of the galaxy after the products from both 
SNe II and SNe la were well mixed. 
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• We found out that the O and Fe radial mass distributions are similar between NGC 4636 
and NGC 1399. This indicates that both objects have similar metal diffusion processes. 
However the MLR of NGC 4636 is 2~3 times larger than that of NGC 1399, possibly due 
to NGC 4636 confining more metal-rich gas. The MLR of NGC 4636 is ~5 times smaller 
than that of groups and clusters of galaxies and there may be a correlation between the 
temperature and MLR in other spherical groups. 

• We confirmed resonance scattering in the Fe xvn line in the central region of NGC 4636, 
first suggested by the XMM-Newton RGS results. The optical depth in the central region 
is, 24.1 < t < 56.3, and thus it is extremely optically thick in the central region if there is 
no turbulence. 

We thank to the anonymous referee for careful reading and many helpful comments. We 
also thank to the Suzaku team for development of hardware/software and operation. 
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